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Abstract. We present a study of the acceleration of electrons at a perpendicular shock that propa- 
gates through a turbulent magnetic field. The energization process of electrons is investigated by uti- 
lizing a combination of hybrid (kinetic ions and fluid electron) simulations and test-particle electron 
simulations. In this method, the motions of the test-particle electrons are numerically integrated in 
the time-dependent electric and magnetic fields generated by two-dimensional hybrid simulations. 
We show that large-scale magnetic fluctuations effect electrons in a number of ways and lead to 
efficient and rapid energization at the shock front. Since the electrons mainly follow along magnetic 
lines of force, the large-scale braiding of field lines in space allows the fast-moving electrons to in- 
teract with the shock front and get accelerated multiple times. Ripples in the shock front occurring 
at various scales will also contribute to the acceleration by mirroring the electrons. Our calculation 
shows that this process favors electron acceleration at perpendicular shocks. The acceleration effi- 
ciency is critically dependent on the turbulence amplitude and coherence length. We also discuss the 
implication of this study for solar energetic particles (SEPs) by comparing the acceleration of elec- 
trons with that of protons. Their correlation indicates that perpendicular shocks play an important 
role in SEP events. 
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INTRODUCTION 

The acceleration of charged particles at collisionless shocks is a remarkable process in 
astrophysics and space physics. The diffusive shock acceleration (DSA) has been very 
successful in explaining the acceleration process and is thought to be the mechanism for 
accelerating solar energetic particles (SEPs). However, how a population of low energy 
particles gain energy, the so-called "injection problem", is still not well understood. 
The acceleration of electrons at collisionless shocks is considered to be more difficult 
than that of ions. This is primarily due to the fact that the gyroradii of electrons are much 
smaller compared with that of protons at the same energy (by a factor of yjnti/m e ~ 43), 
therefore low energy electrons can not resonantly interact with large-scale turbulence or 
ion-scale waves. In contrast, energetic electrons are often observed to be associated with 
quasi-perpendicular shocks. 

In SEP events, energetic electrons are found to be associated with energetic ions. Hag- 
gerty & Roelof [1] reported that energetic electrons (175-315 keV) are well associated 
with protons with energy between 1.8 and 4.7 MeV. Cliver [2] showed that, the electrons 
(~ 0.5 MeV) and the protons (> 10 MeV) in large SEP events are strongly correlated. 
The correlation still exists even when the isotope ratio (Fe/O) changes. The tight cor- 



relation between electrons and ions has also been found in ground base level events 
(Tylka, private communication [3]), where both electrons and protons are accelerated to 
relativistic energies. It is interesting to note that the high energy electrons and protons 
are also closely correlated in solar flares [4]. These observations indicate that a common 
mechanism may exist in the process of particle acceleration during SEP events. 

Particle acceleration at perpendicular shocks has been considered as a fast and effi- 
cient process [5, 6]. Giacalone [7, 8] has performed numerical simulations for the ac- 
celeration of protons at fast shocks that propagates through a plasma that contains large 
scale magnetic fluctuations. The acceleration to high energy is very efficient, which indi- 
cates that there is no injection problem. Recently, Jokipii & Giacalone [9] have proposed 
a mechanism to solve the injection problem that does not require strong pitch-angle 
scattering from small-scale fluctuations. The fast-moving electrons can move along me- 
andering magnetic field lines and travel back and forth across the shock front, and, 
therefore gain energy from the difference between the upstream and downstream flow 
velocities. Using self-consistent hybrid simulations combined with test-particle simula- 
tions for electrons, Guo & Giacalone [10] have found efficient electron acceleration at 
perpendicular shocks moving through a plasma containing large-scale pre-existing up- 
stream magnetic turbulence. The turbulent magnetic field leads to field-line meandering 
that allows the electrons to get accelerated at the shock front multiple times, like in the 
Jokipii & Giacalone picture. In a more recent paper, Guo & Giacalone [11] demonstrated 
that in solar flare region, a low Mach number shock can also efficiently accelerate both 
electrons and ions. In this paper, we present the results of the acceleration of electrons at 
a perpendicular shock that propagates through a turbulent magnetic field. The energiza- 
tion process of electrons is investigated by utilizing a combination of hybrid simulations 
and test-particle electron simulations. We also discuss the implication of this study for 
SEPs by comparing the acceleration of electrons with that of protons. Their correlation 
indicates that perpendicular shocks play an important role in SEP events. 

NUMERICAL METHOD 

We implement a combination of a 2-D hybrid simulation to model the fields and plasma 
flow and a test particle simulation to follow the orbits of electrons. In the first step, we 
employ a two-dimensional hybrid simulation that includes pre-existing large-scale tur- 
bulence [8]. In the hybrid simulation model, the ions are treated fully kinetically and 
thermal (i.e., non-energetic) electrons are treated as a massless fluid. This approach is 
well suited to resolve ion-scale plasma physics that is critical to describe supercritical 
collisionless shocks. We consider a two-dimensional Cartesian grid in the x — z plane. 
All the physical vector quantities have components in three directions, but depend spa- 
tially only on the two variables. A shock is produced by continuously injecting plasma 
from one end (jc = 0) of the simulation box, and reflected elastically at the other end 
(x = L x ). This boundary is also assumed to be a perfectly conducting barrier. The pileup 
of density and magnetic field creates a shock propagating in the — x direction. To include 
the large-scale magnetic fluctuations, a random magnetic field is superposed on a mean 
field at the beginning of the simulation and is also injected continuously at the x = 
boundary during the simulation. The simplified one-dimensional fluctuations have the 



form B(z,/) = 8B(z,t) + Bi, where Bi is the averaged upstream magnetic field. The 
fluctuating component contains an equal mixture of right- and left-hand circularly po- 
larized, forward and backward parallel-propagating plane Alfven waves. The amplitude 
of the fluctuations is determined from a Kolmogorov-like power spectrum: 

P(k) 1 
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in which L c is the coherence scale of the fluctuations. For the simulations presented in 
this study, we take L c = L z , where L z is the size of simulation box in the z direction. As 
a standard case of this study, we consider a turbulence variance o = 8B 2 /B\ = 0.3. The 
size of the simulation box for most of situations is L x x L z = 400c/ (O p i x \02Ac/(O p i, 
where c/ (o p t is the ion inertial length. The Mach number of the flow in the simulation 
frame is Mao = 4.0, and the averaged Mach number in the shock frame is about 5.6. The 
averaged shock normal angle (0b„) = 90°, but we also present the cases for (0b„) = 60° 
and 75°. We also discuss the effect of different values of turbulence variances and 
correlation lengths. The other important simulation parameters include electron and ion 
plasma beta /3 e = 0.5 and /3, = 0.5, respectively, grid sizes Ax = Az = 0.5c/ (O p u time step 
At = O.OIQ^ 1 , the ratio between light speed and upstream Alfven speed c/vai = 8696.0, 
and the anomalous resistivity is taken to be r\ = 1 x \§~ 5 4%(D~1 . The initial spatially 
uniform thermal ion distribution was generated using 40 particles per cell. Under these 
parameters, the upstream Alfven speed is about 34.5 km/s and the shock speed is about 
193 km/s in upstream frame. 

In the second part of our calculation we integrate the motion equation of an ensem- 
ble of test-particle electrons in the electric and magnetic fields obtained in the hybrid 
simulations. We assume non-relativistic motions, which is reasonable because the high- 
est energy of electrons in our study is still non-relativistic. We release a shell distri- 
bution of electrons with energy of 100 eV, which corresponds to an electron velocity 
V e = 30.1U\ = 5Jv t he in the upstream frame, where U\ is upstream bulk velocity in the 
shock frame and v t h e is the thermal velocity of fluid electrons considered in the hybrid 
simulations, respectively. We use a second-order spatial interpolation and a linear tempo- 
ral interpolation, which ensure the smooth variations of the electromagnetic fields. The 
test-particle electrons are released uniformly upstream when the shock has fully formed 
and is far from the boundaries. The numerical technique used to integrate electron trajec- 
tories is the so-called Bulirsh-Stoer method. It is highly accurate and conserves energy 
well. It is fast when fields are smooth compared with the electron gyroradius. The algo- 
rithm uses an adjustable time-step method based on the evaluation of the local truncation 
error. The time step is allowed to vary between 5 x 10~ 4 and 0.1Q.~ e l , where Q. ce is the 
electron gyrofrequency. The ratio £l ce /£l c i is taken to be the realistic value 1836. The 
total number of electrons in the simulation is 1.6 x 10 6 . The electrons that reach the left 
or right boundary are assumed to escape from the shock region and are removed from 
the simulation. The boundary condition in the z direction is taken to be periodic. The 
readers are referred to [12, 10] for more details on the numerical methods. 




zoomed in on the right bottom and a profile is illustrated on the right upper panel. The shock surface is 
shown to be rippled and irregular in different scales. 

SIMULATION RESULTS 

Figure 1 shows a snapshot of the z component of the magnetic field, B z /B\, at t = 
\\0Q. d l . At this time, the shock is fully developed. The position of the shock front is 
clearly seen from the boundary of the magnetic field jump. On the right bottom a small 
region of the simulation domain is zoomed in, which shows small scale irregularities at 
the shock front. Because of the effect of the large-scale turbulence, the shock surface 
become irregular on a variety of spatial scales from small-scale ripples to large-scale 
structures caused by the interaction between the shock and the upstream turbulence 
[7, 13, 10]. Locally, the structure of the shock, shown in top right panel, is still clearly 
a quasi-perpendicular shock. The meandering of field lines close to the shock surface 
helps to trap the electrons at the shock, leading to efficient acceleration. The shock 
ripples also contribute to the acceleration by mirroring electrons between them. This 
is shown by examining the trajectories of some electrons as illustrated in Figure 2. In 
this figure, the top left plot displays the trajectory of a representative electron in the x — z 
plane, overlapped with a 2-D gray-scale representation of B z at Q, c it = 89.0. The upper 
right plot shows the position of this electron (in x) as a function of time. The electron 
bounces back and forth between the ripples for several times. The energy change as a 
function of the x position, corresponding to these reflections is shown in the bottom left 
panel. We find that there are jumps in energy at each of the reflections. The panel on 
the bottom right shows the electron energy as a function of time that also illustrates the 
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FIGURE 2. An electron trajectory analysis. The top left panel shows the trajectory of the electron in 
the x — z plane, overlapped with a contour of B z magnetic field; The top right panel shows the position of 
the electron in the x direction as a function of time; The bottom left panel illustrates the electron energy 
E/Eq as a function of x\ The bottom right panel shows the dependence of electron energy E/Eq on time. 



features of multiple accelerations related to the multiple reflections. Since the electrons 
mainly follow along the magnetic lines of force, the large-scale braiding of field lines in 
space allows the fast-moving electrons to interact with the shock front multiple times. 
The electrons can also travel along the meandering magnetic field and gain energy from 
the velocity difference between upstream and downstream [9]. 

We now discuss the effect of varying the averaged shock normal angle, the value of 
turbulence variances, and the values of correlation lengths. In the top panels of Figure 3 
we present the electron energy spectra for left: three different mean shock-normal angles 
((Obh) = 60°, 75° and 90°, respectively), middle: three cases with different turbulence 
variances o = 0. 1, 0.3, and 0.5, respectively, and right: three cases with different sizes of 
the simulation box and correlation length L z = 2000, 1024, and 400c/ CO p {. All the spectra 
are obtained in the end of the simulation Q. c jt = 120 and do not change any more. We 
find this process prefer a perpendicular shock, in which more electrons are accelerated 
to high energy. It is also found that the energy spectra harden at high energies when the 
turbulence variance and/or turbulence correlation length is larger, which indicates that 
the large-scale turbulence is important for accelerating electrons to high energies. The 
reason is that the meandering of field lines is enhanced, which allows the electrons have 
a better chance to travel though the shock multiple times. 

We explicitly compare the energy spectra of electrons with that of protons. In the 
bottom panels of Figure 3 we show the energy spectra of downstream protons from the 
hybrid simulations, corresponding to the spectra of electrons in the top panels. We find 




FIGURE 3. Top: energy spectra of accelerated electrons from test-particle electron simulations. Bot- 
tom: The energy spectra of protons directly from hybrid simulations. See description in the text for details. 



more efficient acceleration for protons can be obtained in the case of larger values of 
turbulence variances and correlation lengths. These agree well with the characteristics 
of the acceleration of electrons. In our simulations, the accelerated protons at the oblique 
shock with = 60° are found to reach higher energies, which is different than 

previous works [5, 6, 7, 8]. This is probably due to the limited temporal and spatial scales 
of our simulations. As shown by Giacalone [8] using test-particle simulations, the energy 
spectra of protons reach the highest energy in perpendicular shock case in a longer time 
scale Q. c jt ~ 50000 (this corresponds to 5 — 10 minutes for typical parameters in solar 
corona). However the current results from hybrid simulations do show a population 
of thermal protons can be accelerated to high energies in perpendicular shocks, which 
supports the idea that both electrons and protons can be efficiently accelerated by shocks 
with large shock normal angles. 



SUMMARY 

We have presented the results of the acceleration of electrons (and also protons) at a 
perpendicular shock that propagates through a turbulent magnetic field. The acceleration 
of electrons are enhanced due to the effect of large-scale magnetic turbulence. The 
accompanying results for protons qualitatively show the correlation between accelerated 
electrons and accelerated ions in oblique shocks with large shock normal angles. Since 
the tight correlation between electrons and ions has been observed in SEP events. This 



result indicates that quasi-perpendicular/perpendicular shocks play an important role in 
SEP events. 
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